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Anisotropic effect on dynamics of block copolymers in lamellar 
phases: Relaxation and the grain boundary motion 
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We consider the effects of anisotropic diffusion and hydrodynamic flows on the relaxation time scales of the lamellar phase of a 
diblock copolymer. We first extend the two-fluid model of a polymer solution to a block copolymer, and include a tensor mobility 
for the diffusive relaxation of monomer composition which is consistent with the uniaxial symmetry of the lamellar phase. The 
resulting equation is coupled to the momentum conservation equation, allowing also for a dissipative stress tensor for a uniaxial 
fluid. We then study the linear relaxation of weakly perturbed lamellae, and the motion of a tilt grain boundary separating two 
semi-infinite domains. We find that anisotropic diffusion has a negligible effect on the linear relaxation of the layered phase (in 
the long wavelenght limit), whereas the introduction of hydrodynamic flows considerably speeds the decay to a rate proportional 
to Q 2 , where Q <C 1 is the wavenumber of a transverse perturbation to the lamellar phase (diffusive relaxation scales as Q 4 
instead). On the other hand, grain boundary motion is siginificantly affected by anisotropic diffusion because of the coupling 
between undulation and permeation diffusive modes within the grain boundary region. 



1 Introduction 

Block copolymers are finding numerous applications in nan- 
otechnology 1 , and have been of great interest in soft matter 
science because they undergo microphase separation to or- 
dered phases of different symmetries (For a brief review see 
Ref. 2 and references therein). Order parameter models to de- 
scribe equilibrium properties and the microphase phase di- 
agram were given in Refs.- - -. Microphase separation also 
brings about interesting dynamical properties, including un- 
usual rheological response^ - - and orientation selection during 
shear aligning .J£JJ- 

■ In principle, the low-frequency and long-wavelength hydro- 
dynamic equations of motion for ordered systems can be de- 
— H .rived by considering conservation laws and symmetry argu- 
^ ments.i 2 ^ For a diblock copolymer there are six conserva- 
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tion laws: Total mass of each monomer, three components 
of momentum, and energy. The two-fluid model has been 
widely used to describe the dynamical behavior of polymer 
solutions and blends in this hydrodynamic regimeJ 4 -"— . More 
recently, the two-fluid model has been further extended to de- 
scribe the dynamics of block copolymer melt s 17 i 18 . The result- 
ing governing equations involve dissipative dynamics for the 
order parameter that represents local monomer composition 
changes, and overdamped (or Stokesian) dynamics for the ve- 
locity fields. However, both dissipative relaxation of the order 
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parameter and hydrodynamic equations assume that the block 
copolymer phase is an isotropic fluid. For example, the order 
parameter equation contains a current associated with free en- 
ergy dissipation due to the relative motion of the two types of 
monomers, with a kinetic coefficient (or mobility) that does 
not respect the uniaxial symmetry of the lamellar phase. 

We address here the constitutive relations for a block 
copolymer lamellar phase that reflect the uniaxial symmetry 
of this ordered phase, and derive the corresponding equations 
of motion in the hydrodynamic regime (Sec. II). This requires 
both an anisotropic mobility tensor for order parameter diffu- 
sion, and anisotropic viscosities for the dissipative stress ten- 
sor. We next investigate in Sec. Ill the linear relaxation of 
weakly perturbed lamellae, and in Sec. IV the motion of a tilt 
grain boundary. 

2 Model equations 

The equations governing the dynamics of an AB -diblock 
copolymer in its lamellar phase can be derived by using con- 
servation laws and broken symmetry argument s 12 i 1 3 i 19 . Ac- 
cording to the two-fluid model for polymer solutions, blends, 
or diblock copolymers^ 4 - - — one introduces two continuity 
equations for each monomer 

dtp A = -di(p A ViA), (1) 
dtpB = -di(pBViB), (2) 

where pa, b and v^, b are the monomer number fraction and 
the corresponding velocities with A and B denoting polymers 



This journal is ©The Royal Society of Chemistry [year] 



Journal Name, 20 1 0, [vol], 1 -0 | 1 



and solvent for polymer solutions or two types of monomers 
in polymer blends or diblock copolymers. It is customary to 
introduce an order parameter i/j = pA — Pb so that 



dti) + diJi = 0, 



(3) 



where the flux Ji has both a reversible part Jf that accounts 
for advection of ip, and a dissipative part Jf responsible for 
energy dissipated due to the relative motion of the two types 
of monomers. By taking the difference of the continuity equa- 
tions, the order parameter equation Eq. ® contains a re- 
versible current Jf = tpVi where the average flow velocity 
is Vi = phViA + PbVib with pa + Pb = 1, and the dissipative 
current 



jf = 2pApB 



ViA ~ V iB 



(4) 



According to the two-fluid model, the relative velocity be- 
tween two monomers in Jf is obtained by employing a 
Rayleigh's variational principle 20 in which a dissipative func- 
tion is introduced that contains the square of the relative ve- 
locity times an isotropic friction coefficient (. Under this as- 
sumption, the resulting dissipative current is-^~— 



Jf = -cWb) 2 ^, 



(5) 



where p = pa — Mb is the difference in the monomer chemi- 
cal potentials. The constitutive law, however, must reflect the 
ordered phase's symmetry. In general, one would write 

Jf = -Aijdjfj,, 



with Aij an anisotropic kinetic constant tensor. The number 
of independent components of is determined by the sym- 
metry of the phase. For layered systems of uniaxial symmetry, 
there are only two independent components, 



A, 



A L nirij + A T (Sij - riirij), 



(7) 



where ni is the unit normal to the layers. In the hydrodynamic 
limit, and when lamellae are weakly perturbed, the unit vec- 
tor rii varies slowly. We assume that in this limit A^ retains 
the same form as above with the local principal axis defined 
by the local normal to the disturbed layers. Furthermore, one 
would expect, and this is confirmed experimentally, that since 
chain mobility across the layers is suppressed relative to mo- 
tion parallel to the layers because of large entropic barriers for 
copolymers to move across layers,— then Al < At. 

For lamellar diblock copolymers in the weak segregation 
limit the relative chemical potential is obtained by taking the 
functional derivative with respect to ip of the free energy func- 
tional given by-~- 



d 6 r< 



V- 



(V 2 + g 2 )^ 



>, (8) 



where qo is the wavenumber of the layers, and K, and £ are 
coefficients that depend on the material properties of the block 
copolymer. Since the reversible current Jf = ipVi, and with 
Eq. © for Jf, the order parameter equation Eq. © becomes 



d t ip + di (ijjVi) = di 



Aijdj 



6£ 

5i/j 



(9) 



Given the high viscosity of block copolymer melts, the 
momentum conservation equation is considered in the over- 
damped limit (small Reynold number) 19 



d z P 



d 3 af 3 = 0, 



(10) 



which has an implicit dependence on the velocity. In this equa- 
tion P is the pressure, afj the reversible elastic stress tensor, 
and afj the dissipative stress tensor. The gradient of the re- 
versible elastic stress tensor is simpl y 19 i 22 i 23 
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For systems with uniaxial symmetry there are five indepen- 
dent viscosities in <jP- 24 , although the number reduces to three 
under the assumption of incompressibility. Hence the dissipa- 
tive stress tensor for an incompressible system with uniaxial 
symmetry can be written as^i 

ofj = a^injnkniVki^a^Vij^a^nkiniVkj^rijVki), (12) 

(6) 

where the strain rate tensor is Vij = (d{Vj + djVi)/2. Again, 



we will assume that for slowly varying lamellar phases, the 
expression (fT2l) holds locally. 

Finally, block copolymers are normally assumed to be in- 
compressible fluids di Vi =0. Consequently, we have Eqs. ©- 
(fT2l) plus the incompressibility condition diVi = as the gov- 
erning equations for the evolution the lamellar phase of a di- 
block copolymer. 

Before we proceed any further, we recast the governing 
equations of motion, Eqs. (l9b-([T2b. in terms of dimensionless 
quantities x' = g x, t' = ^A L qQt, t/j' = ip/^/£/xQo> and 
F' = F / (t; 2 q$ / x) • With the newly defined variables the or- 
der parameter equation can be rewritten as 



d t ip + Vidiip di 



with the free energy being 



Aijdj 



5£ 

5ip 



= - I d 3 H 



(V 2 + g 2 )V> 



and the anisotropic kinetic tensor 



A* 



tiitij + A(6ij - riiUj), 



(13) 



(14) 



(15) 
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where e = K/^q 4 and A = At/A^. In addition, the momen- 
tum conservation equation becomes 



0, 



(16) 



where we have used the viscosity scale to rescale the dissi- 
pative stress tensor such that 

of. = ainifijUkniVki + + a 56 n k (niV k j + njV ki ). (17) 

In the rescaled momentum conservation equation a coeffi- 
cient C(= £goM£X a 4) appears because of the rescaling. The 
reader should also note that we have omitted the primes for 
clarity, and retained qo explicitly in the free energy, although 
it becomes unity in the rescaled units. 

3 Relaxation of Diblock Copolymer Lamellae 

In order to ascertain the effects of an anisotropic diffusivity 
and the coupling to hydrodynamics flows, we begin by investi- 
gating the linear relaxation of weakly perturbed lamellae. This 
linear analysis can be done analytically, and provides us with 
relaxation times that depend on the strength of the anisotropy 
and the hydrodynamic coupling coefficient. Let us consider 
as reference state a stationary solution of the order parameter 
equation Eq. (TT3l) for a lamellar phase, 



compared to the amplitude ipi of the reference state. We fol- 
low Zhang's study— of the effect of hydrodynamic flows on 
the relaxation of lamellar block copolymers to derive the am- 
plitude equations by replacing Eq. (l22l) into Eq. (12T1) and by 
retaining only terms linear in ^2,3 and 3 . 

By substituting the perturbation Eq. (E21) into Eq. (I2TK the 
term in squared brackets becomes, in Fourier space (k, oj), 



SF 



(k,w) 

s-(k 2 
d 3 kidooi 



9o 2 ) 2 



V>(k,w) 

d 3 k2du)2 

x ^(k 



^(ki,cji)^(k 2 ,cJ2)x 

- ki — k 2 , (jJ - uj\ - UJ2) 
(2ir) 4 MoS(uj)5(k + q) + (27r) 4 M 5(uj)S(k - q) 



(2-) 2 



Mi(o;)(5(k - q - Q) + MZ(u)5(k + q + Q) 



• M 2 (u)5(k - q + Q) + M*(u)5(k + q - Q) 



where 



M 



1>i 



e-(q'-ql?--^l 



(23) 
(24) 



^ s (r) = ^1 cos(q- r). 



(18) 



The wave vector q defines the unit normal to lamellae, rii = 
qi I q, and the kinetic constant tensor associated with the refer- 
ence state is 



Mi(cj) 



e^ 2 {uj) + 
3 



q + Q| 2 -9o 

2^H + V 3 *M 



(25) 



(19) 



Since A^ is uniform for the reference lamellae ^ s , the order 
parameter equation, Eq. ([T3l) , reduces to 



d t ip + Vidiip = Aijdidj ( — 



fSF 



or 



(20) 



. (21) 



We now consider small disturbances of wave number (Q <C (7) 
such that 



^(r, t) = ^1 cos(q • r) + ip 2 (t) exp[i(q + Q) • r] 
+ fait) exp[z(q - Q) • r] + c.c. + . . . , 



(22) 



where c.c. stands for complex conjugation, and the perturba- 
tive amplitudes ip 2 , ^3 an d their complex conjugates are small 



M 2 (u) 



1 2 



|q-QI : 



2^ 3 (w) + ^(w) 



(26) 



where 3 (u) are the temporal Fourier transforms of the com- 
plex conjugates of 7/>2,3 (t), respectively. Then we obtain the 
R.H.S of Eq. (ED by multiplying Eq. (EU) by -Kijhkj. 

Since there is no external source to generate hydrodynamic 
flows, the flow velocity v is produced by changes in ip, and 
solely determined by the momentum conservation equation. 
In order to obtain the velocity in terms of ip we take the gra- 
dient of the momentum conservation equation Eq. dT6b . and 
solve for the pressure in Fourier space. Next, the obtained 
pressure is substituted back into the momentum conservation 
equation, resulting in 



2C 



1 



1 + a 56 k 2 



Sij 



ki kj 



(k,w). 



(27) 
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It is worth mentioning here that the flow velocity v\ does not 
depend on the viscosity coefficient ol\ because of the incom- 
pressibility condition. We now obtain the velocity in terms 
of the amplitudes of the perturbation by substituting Eqs. (l22l) 
and ([23]) into Eq. (l27lh 



= -i(27r) 3 



2C q 3 



1 



Sin 



^ ^(k + Q)x 



®56 k 2 \ k 2 



2 

i(27r) 3 



2C 



2 



a 56 k 2 



Sij 



(28) 



in which we have retained only ±Q modes because the flow 
velocity is already linear in the perturbation amplitudes, and in 
the convective term of the order parameter equation it couples 
only to the reference state ip s of modes ±q. Thus we find for 
the convective term in the order parameter equation, 



/ 



(Pkxdujx 

Vi(k - ki,o; - cji)zfeii^(ki,a;) 

(2tt) 4 



^(k + q,cj) - Ui(k- q,uj) 



(29) 



Finally combining Eqs. (l23t and (l29l) we obtain the am- 
plitude equations, and present them separately for different 
modes in the small Q limit. First there is an equation, cor- 
responding to modes ±Q, which defines the amplitude of the 
reference wave ipi 



2\2 



(Q 2 ~ Q 2 o) 



(30) 



It is required that e > (q 2 — q^) 2 to ensure that ^ s exists. 
Second, there are two amplitude equations corresponding to 
the modes q + Q and — q + Q for ^2 and ^3 



a b 
c d 



where 
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(31) 

(32) 
(33) 
(34) 



3?/; 2 

d= -iw + (JJ + A_)Z_ + ^-A_. (35) 
We have introduced a hydrodynamic coupling coefficient 



H= 2 - 



c 



1 



3 1 + a 56 Q 2 
and have defined 



qiqj 



(tf 1 - & 



A ± =q 2 ±2(q l Q l ) + ^^+A 



Q 2 



i± = (|q±Q| a -(zg) a -((Z a -(zg) a . 



(36) 



, (37) 



(38) 



Additionally, there are two amplitude equations for modes 
— q — Q and q — Q leading to equatoins for for ^| and ips 
that are the complex conjugates of Eq. (|3TT) . 

In general, the perturbation is linearly stable if the fre- 
quency uj is negative and pure imaginary. The imaginary fre- 
quencies are the decay rates, and they are obtained by solving 
the characteristic equation of the 2 x 2 matrix in Eq. (I3T1) . Here 
we consider two different cases: an undulation mode (Q _L q) 
and a permeation mode (Q || q). When Q || q, from Eqs. (l36t 
and (l37t we find that neither anisotropic diffusion nor hydro- 
dynamic coupling affect the relaxation rate. From Eq. (|3TT) we 
find two relaxation rates 



32, 



+ 2(11^ - 9q£)q 2 + -{q 2 - q^ 
0(Q% 



Q 2 



(39) 



4oV-f(4 2 



Q 2 



(40) 



+ 0(Q% 

The relaxation rate describes the decay of the perturbation 

amplitude, whereas describes decay of its phase. Due to 

the order one term in the amplitude decays much faster 
than the phase, and follows adiabatically any change in the 
phase. 

When Q _L q, we have l± = 2(q 2 - ql)Q 2 + Q 4 , and 
A± = q 2 + AQ 2 . In this case Eq. (PIT) reduces to 



a_L 



b± 



0, 



where 



(# + A + )Z + + ^A + , 



-S7+ 



(41) 

(42) 
(43) 
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Then it is straightforward to obtain two relaxation rates 

3 



2(^ 



Q 2 + 0(Q 4 ), (44) 



rl\ =4h(q 2 - ql) + 2 



h + (q Z - ql)q 2 



Q 2 



q 2 + 2A(q 2 - ql) 



Q 4 + 0(Q 6 ), 



where 



c 



3 1 + a 56 



q 



(Q 2 - <ll? 



(45) 



(46) 



Again these two rates rj\ and rj \ govern the relaxation of 
the amplitude and phase, respectively. We find that anisotropic 
diffusion contributes at order Q 2 to rj\, and at order Q 4 to 
rj 2 '•> therefore, its effect on the relaxation of weakly perturbed 
lamellae is negligible in the limit of small Q. Hydrodynamic 
flow only couples to the phase of perturbation, and derives the 
phase to decay faster for q ^ q unlike the previous case of 
parallel perturbations. If q = go, from Eq. (l46l) we can define 
the hydrodynamic diffusion length Ahydro as 



\2 

^hydro 



1 

hql 



3 A L (q 4 
2 



«56)X 



(47) 



in the dimesional units. Then, the marginal mode for instabil- 
ity, the diffusion due to the hydrodynamic flow becomes neg- 
ligible with respect to the isotropic order parameter diffusion 
when Ahydro Q > V%- In this case the diffusion mode becomes 
identical to the undulation mode of Smectic-A liquid crystal 
(~ Q 4 )- 25 When AhydroQ < V%, the hydrodynamic diffusion 
dominates with a decay rate proportional to Q 2 . 

4 Grain Boundary Motion 

We have studied in the previous section the effect that an 
anisotropic diffusion tensor can have on the linear relaxation 
dynamics of a weakly perturned lamellar phase. In general, 
however, when a block copolymer is brought below its mi- 
crophase separation transition point, a large number of struc- 
tural defects such as grain boundaries are quenched in a spa- 
tially extended system. These boundaries separate locally lay- 
ered domains of different orientations producing a macroscop- 
ically inhomogeneous system. Within the defected region (of 
extent which is much larger that the lamellar wavelength in 
the weak segregation limit considered in this paper) undula- 
tion and permeation diffusive modes strongly couple. This 
coupling leads to a significant contribution to defect dynam- 
ics, the subject matter of this section. 

Grain boundary motion in layered systems has been inten- 
sively studied in other contexts such as Rayleigh-Benard con- 
vection^"—. In the governing equation for Rayleigh-Benard 



Region III Region II Region I 




Fig. 1 Sketch of a 90° tilt grain boundary separating two 
semi-infinte lamellae. In regions I and III the lammellae are nearly 
in their stationary state with wavevectors qox and qoy, respectively. 
The grain boundary region (II) is the region between the two blue 
dashed lines. In it, both envelopes A and B exhibit large variation. 



convective rolls, the order parameter is the vertical veloc- 
ity which is equivalent to the order parameter in equation 
Eq. (Tl3l) . In this section we examine the consequences of al- 
lowing an anisotropic kinetic coefficient, compatible with the 
uniaxial symmetry of a layered block copolymers. Our results 
extend the analysis of Rayleigh-Benard convective rolls given 
in Refs.^2. 

We consider a 90° tilt grain boundary that separates 
two semi-infinite domains of block copolymer lamellae with 
wavevector g x m region I, and qoy in region III as shown 
in Fig. [T] For simplicity, we focus on an effective two- 
dimensional system by taking advantage of translational sym- 
metry of lamellae along the direction (z), perpendicular to the 
wavevectors of the two semi-infinite lamellae. 

In order to take into account the inhomogeneous nature of 
the kinetic constant tensor due to the presence of a 90° 
grain boundary in the interface region II, we model with an 
auxiliary function Q(x) that interpolates smoothly within the 
width of the grain boundary from unity in the region where one 
semi-infinite lamella is present to zero in the opposite region: 



6(x)A 1 i7 +e(-x)Al J 



@{x) 



ASij + (1 - A)^^ 



+e(-a?) 



A^- + (1 



A)S iy S jy 



(48) 



where A^ and A^ 1 are the kinetic constant tensors defined in 
regions I and III with q$x and q^y, respectively, and Q(—x) = 
1 — 0(x). In writing A^ as Eq. (|48t we have assumed that 
the two semi-infinite lamellae decay with a common length 
scale proportional to the size of grain boundary. In Refs.^2i 
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it is found that the width of grain boundary diverges as e -1 / 2 
in the weak segregation limit, so that we can infer that 
changes very slowly in the interfacial region. 

Due to the inhomogeneity of its gradient does not van- 
ish, and the order parameter equation becomes 



(49) 



where the coupling to hydrodynamic flow is neglected for the 
analysis in this section. We now use a multiple scale analy- 
sis to derive amplitude equations close to the linear instability 
threshold (e <C 1). Following Tesauro and Cross 27 , we intro- 
duce slow variables X = e 1 / 4 ^, X = e l / 2 x, Y = e 1 ^ A y, 
Y = e 1 / 2 ?/, T = et, and expand the derivatives d x —> 
d x + e^ 4 d x + e^ 2 d x , d y -> d y + e l ^d Y + e^dy, and 
d t —> edr- Next, by noting that the order parameter scales as 
e 1 / 2 we take 



V3 



Aex.-p(iqox) + B ex.p(iqoy) 



c.c, 



(50) 



where A and B are functions of the slow variables only. Then 
the amplitude equations at C?(e 3 / 2 ) are, 



-it 



SF' 
<L4* 



iqod x O(x) + q%®(-x) 



SF' 
SA*' 



(51) 
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d T B = - q 2 — + Aq 2 e(x)^, 



where the free energy functional is 

F'[A,A*,B,B*] 



/ 



d 2 rl - \A\ 2 - \B\ 2 + x - ( |A| 4 + |5| 4 ) + 2\A\ 2 \B\ 2 



(2iq d x + d 2 )A 



(2iq d y + d 2 x )B 



(53) 



and A = 1— A. Since A is positive, it is required that A < 1. It 
is easy to show that the isotropic case of dissipative coefficient 
(Aij ~ Sij) is recovered by taking A = 1. 

The solutions A s and B s of the stationary planar 90° grain 
boundary without anisotropy, given in Refs.— and—, are the 
stationary solutions of the amplitude equations. Both A s and 
B s saturate to e 1 / 2 as x tends to -f oo or — oo respectively, but 
have different decaying behaviors within the grain boundary. 
The amplitude A s has a longer decaying length scale ~ e -1 / 2 
than B s . Since the width of the grain boundary scales as e -1 / 2 , 



it is reasonable to assume that this is the same scale of vari- 
ation of the function Q(x). This leads to c^A^- ~ e 1 / 2 , and 
its contribution to the amplitude equations appears at higher 
order in e. The remaining relevant term is Aijdidj(5F/5ip), 
and Eqs. (TSTb and (f52l) reduce to 



d T A 



d T B 



1 - AB(-.t) 



AQ(x) 



SF' 

SF' 
SB*' 



(54) 



(55) 



We now use the energy method^* 2 ^ to calculate the grain 
boundary velocity. The time derivative of F' is given by 



dF' 
~dt 



d 2 



-\dtAf 



q S j 1 - A9(-x) 1 
q 2 J 1 - AB(x) 1 1 



(56) 



Note that the anisotropy effect is on the R.H.S. only, and that 
the coefficients of both |c^A| 2 and \dtB\ 2 are positive because 
A < 1. When the grain boundary moves with a velocity ^gb, 
it is convenient to take A s (x — vbg^) and B s (x — ^bg^) so 
that they are stationary in the moving frame. Then the time 
derivative can be replaced d t with —vbg9 x , and we find for 
the grain boundary velocity 



v GB = M / dy 



T'{x = oo, y) — T'{x = —oo, y) 



(57) 



where T' is the free energy density, and the effective mobility 
of the boundary is given by 



1 

M 



d 2 r 



1 



-\d x A 



s|2 



</o 2 



l-Ae(-i) 1 



Since A is an order one quantity, the contribution to the 
boundary velocity due to anisotropic diffusion is large. The 
envelope A relaxes in region I (of dominant orientation q$x) 
differently than in region III (of dominant orientation qoy). 
In region I, diffusion is along the lamellar normal, whereas 
this component of the order parameter evolves through trans- 
verse difussion in region III. Exactly the same is true of com- 
ponent B. As a consequence, the boundary velocity depends 
on a weighted average of the two independent diffusion co- 
efficients, with the weight function being the gradient of the 
order parameter envelopes, as given in Eq. (l58l) . Of course, a 
similar qualitative behavior can be expected in the vicinity of 
other structural defects. When A = 0, the isotropic result of 
Refs.^ 2 ^ is recovered. 

We note that there is no grain boundary motion for unper- 
turbed lamellae when T'{x = oo) — T'{x — — oo) ~ e. In 
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practice, an imbalance of the free energies caused by external 
sources is necessary to drive grain boundary motion so as to 
reduce excess free energy. We find that anisotropy enhances 
(reduces) v G b when A<0(0<A<1). 

5 Conclusion 

We have investigated diffusive relaxation in lamellar phases 
of block copolymers when allowing for uniaxial symmetry of 
the consitutive law between diffusive forces and fluxes, as well 
as hydrodynamic coupling. We have shown that coupling to 
flows leads to a relaxation rate proportional to Q 2 , where Q 
is the wavenumber of the characteristic perturbation. The uni- 
axial symmetry of the lamellar phase of a diblock copolymers 
requires an anisotropic kinetic constant in the order parameter 
equation, and an anisotropic stress tensor in the momentum 
conservation equation. With them, we have calculated the re- 
laxation rates of a weakly perturbed lamella, and found that 
the effect of anisotropy becomes negligible compared to ei- 
ther hydrodynamic flow or (isotropic) order parameter diffu- 
sion. We have also studied the motion of a grain boundary by 
calculating its velocity, and shown that the velocity is signifi- 
cantly affected by anisotropic diffusion. 
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